Search for: All records

Abstract This paper conducts a multi‐instrument and data assimilation analysis of the three‐dimensional ionospheric electron density responses to the total solar eclipse on 08 April 2024. The altitude‐resolved electron density variations over the continental US and adjacent regions are analyzed using the Millstone Hill incoherent scatter radar data, ionosonde observations, Swarm in situ measurements, and a novel TEC‐based ionospheric data assimilation system (TIDAS) with SAMI3 model as the background. The principal findings are summarized as follows: (a) The ionospheric hmF2 exhibited a slight enhancement in the initial phase of the eclipse, followed by a distinct reduction of 20–30 km in the recovery phase of the eclipse. The hmF2 in the umbra region showed a post‐eclipse fluctuation, characterized by wavelike perturbations of 10–25 km in magnitude and a period of 30 min. (b) There was a substantial reduction in ionospheric electron density of 20%–50% during the eclipse, with the maximum depletion observed in the F‐region around 200–250 km. The ionospheric electron density variation exhibited a significant altitude‐dependent feature, wherein the response time gradually delayed with increasing altitude. (c) The bottomside ionospheric electron density displayed an immediate reduction after local eclipse began, reaching maximum depletion 5–10 min after the maximum obscuration. In contrast, the topside ionospheric electron density showed a significantly delayed response, with maximum depletion occurring 1–2.5 hr after the peak obscuration. 

Abstract This study investigates the ionospheric total electron content (TEC) responses in the 2‐D spatial domain and electron density variations in the 3‐D spatial domain during the annular solar eclipse on 14 October 2023, using ground‐based Global Navigation Satellite System (GNSS) observations, a novel TEC‐based ionospheric data assimilation system (TIDAS), ionosonde measurements, and satellite in situ data. The main results are summarized as follows: (a) The 2‐D TEC responses exhibited distinct latitudinal differences. The mid‐latitude ionosphere exhibited a more substantial TEC decrease of 25%–40% along with an extended recovery time of 3–4 hr. In contrast, the equatorial and low‐latitude ionosphere experienced a smaller TEC reduction of 10%–25% and a faster recovery time of 20–50 min. The minimal eclipse effect was observed near the northern equatorial ionization anomaly crest region. (b) The ionospheric electron density variations during the eclipse were effectively reconstructed by TIDAS data assimilation in the 3‐D domain, providing important altitude information with validity. (c) The ionospheric electron density variations showed a notable altitude‐dependent feature. The eclipse led to a substantial electron density reduction of 30%–50%, with the maximum depletion occurring around the ionospheric F2‐layer peak height (hmF2) of 250–350 km. The post‐eclipse recovery of electron density exhibited a relatively slower pace near the F2‐layer peak height than that at lower and higher altitudes. 

This paper presents a multi-instrument observational analysis of the equatorial plasma bubbles (EPBs) variation over the American sector during a geomagnetically quiet time period of 07–10 December 2019. The day-to-day variability of EPBs and their underlying drivers are investigated through coordinately utilizing the Global-scale Observations of Limb and Disk (GOLD) ultraviolet images, the Ionospheric Connection Explorer (ICON) in-situ and remote sensing data, the global navigation satellite system (GNSS) total electron content (TEC) observations, as well as ionosonde measurements. The main results are as follows: 1) The postsunset EPBs’ intensity exhibited a large day-to-day variation in the same UT intervals, which was fairly noticeable in the evening of December 07, yet considerably suppressed on December 08 and 09, and then dramatically revived and enhanced on December 10. 2) The postsunset linear Rayleigh-Taylor instability growth rate exhibited a different variation pattern. It had a relatively modest peak value on December 07 and 08, yet a larger peak value on December 09 and 10. There was a 2-h time lag of the growth rate peak time in the evening of December 09 from other nights. This analysis did not show an exact one-to-one relationship between the peak growth rate and the observed EPBs intensity. 3) The EPBs’ day-to-day variation has a better agreement with that of traveling ionospheric disturbances and atmospheric gravity waves signatures, which exhibited relatively strong wavelike perturbations preceding/accompanying the observed EPBs on December 07 and 10 yet relatively weak fluctuations on December 08 and 09. These coordinate observations indicate that the initial wavelike seeding perturbations associated with AGWs, together with the catalyzing factor of the instability growth rate, collectively played important roles to modulate the day-to-day variation of EPBs. A strong seeding perturbation could effectively compensate for a moderate strength of Rayleigh-Taylor instability growth rate and therefore their combined effect could facilitate EPB development. Lacking proper seeding perturbations would make it a more inefficient process for the development of EPBs, especially with a delayed peak value of Rayleigh-Taylor instability growth rate. 

Abstract This study has developed a new TEC‐based ionospheric data assimilation system for 3‐D regional ionospheric imaging over the South American sector (TIDAS‐SA) (45°S–15°N, 35°–85°W, and 100–800 km). The TIDAS‐SA data assimilation system utilizes a hybrid Ensemble‐Variational approach to incorporate a diverse set of ionospheric data sources, including dense ground‐based Global Navigation Satellite System (GNSS) line‐of‐sight Total Electron Content (TEC) data, radio occultation data from the Constellation Observing System for Meteorology, Ionosphere, and Climate‐2 (COSMIC‐2), and altimeter TEC data from the JASON‐3 satellite. TIDAS‐SA can produce a reanalyzed three‐dimensional (3‐D) electron density spatial variation with a high time cadence, yielding spatial‐temporal resolution of 1° (latitude) × 1° (longitude) × 20 km (altitude) × 5 min. This allows us to reconstruct and study the 3‐D ionospheric morphology with multi‐scale structures. The performance of the data assimilation system is validated against independent ionosonde and in situ measurements through an experiment for a strong geomagnetic storm event on 03–04 November 2021. The results demonstrate that TIDAS‐SA can provide detailed and altitude‐resolved information that accurately characterizes the storm‐time ionospheric disturbances in vertical and horizontal domains over the equatorial and low‐latitude regions of South America. 

The Tonga volcano eruption at 04:14:45 UT on 2022-01-15 released enormous amounts of energy into the atmosphere, triggering very significant geophysical variations not only in the immediate proximity of the epicenter but also globally across the whole atmosphere. This study provides a global picture of ionospheric disturbances over an extended period for at least 4 days. We find traveling ionospheric disturbances (TIDs) radially outbound and inbound along entire Great-Circle loci at primary speeds of ∼300–350 m/s (depending on the propagation direction) and 500–1,000 km horizontal wavelength for front shocks, going around the globe for three times, passing six times over the continental US in 100 h since the eruption. TIDs following the shock fronts developed for ∼8 h with 10–30 min predominant periods in near- and far- fields. TID global propagation is consistent with the effect of Lamb waves which travel at the speed of sound. Although these oscillations are often confined to the troposphere, Lamb wave energy is known to leak into the thermosphere through channels such as atmospheric resonance at acoustic and gravity wave frequencies, carrying substantial wave amplitudes at high altitudes. Prevailing Lamb waves have been reported in the literature as atmospheric responses to the gigantic Krakatoa eruption in 1883 and other geohazards. This study provides substantial first evidence of their long-duration imprints up in the global ionosphere. This study was enabled by ionospheric measurements from 5,000+ world-wide Global Navigation Satellite System (GNSS) ground receivers, demonstrating the broad implication of the ionosphere measurement as a sensitive detector for atmospheric waves and geophysical disturbances. 

Abstract A new TEC‐based ionospheric data assimilation system (TIDAS) over the continental US and adjacent area (20°–60°N, 60°–130°W, and 100–600 km) has been developed through assimilating heterogeneous ionospheric data, including dense ground‐based Global Navigation Satellite System (GNSS) Total Electron Content (TEC) from 2,000+ receivers, Constellation Observing System for Meteorology, Ionosphere, and Climate radio occultation data, JASON satellite altimeter TEC, and Millstone Hill incoherent scatter radar measurements. A hybrid Ensemble‐Variational scheme is utilized to reconstruct the regional 3‐D electron density distribution: a more realistic and location‐dependent background error covariance matrix is calculated from an ensemble of corrected NeQuick outputs, and a three‐dimensional variational (3DVAR) method is adopted for measurement updates to obtain an optimal state estimation. The spatial‐temporal resolution of the reanalyzed 3‐D electron density product is as high as 1° × 1° in latitude and longitude, 20 km in altitude, and 5 min in universal time, which is sufficient to reproduce ionospheric fine structure and storm‐time disturbances. The accuracy and reliability of data assimilation results are validated using ionosonde and other measurements. TIDAS reanalyzed electron density is able to successfully reconstruct the 3‐D morphology and dynamic evolution of the storm‐enhanced density (SED) plume observed during the St. Patrick's day geomagnetic storm on 17 March 2013 with high fidelity. Using TIDAS, we found that the 3‐D SED plume manifests as a ridge‐like high‐density channel that predominantly occurred between 300 and 500 km during 19:00–21:00 UT for this event, with the F2 region peak height being raised by 40–60 km and peak density enhancement of 30%–50%. 

Abstract This paper investigates the local and global ionospheric responses to the 2022 Tonga volcano eruption, using ground‐based Global Navigation Satellite System total electron content (TEC), Swarm in situ plasma density measurements, the Ionospheric Connection Explorer (ICON) Ion Velocity Meter (IVM) data, and ionosonde measurements. The main results are as follows: (a) A significant local ionospheric hole of more than 10 TECU depletion was observed near the epicenter ∼45 min after the eruption, comprising of several cascading TEC decreases and quasi‐periodic oscillations. Such a deep local plasma hole was also observed by space‐borne in situ measurements, with an estimated horizontal radius of 10–15° and persisted for more than 10 hr in ICON‐IVM ion density profiles until local sunrise. (b) Pronounced post‐volcanic evening equatorial plasma bubbles (EPBs) were continuously observed across the wide Asia‐Oceania area after the arrival of volcano‐induced waves; these caused aN_edecrease of 2–3 orders of magnitude at Swarm/ICON altitude between 450 and 575 km, covered wide longitudinal ranges of more than 140°, and lasted around 12 hr. (c) Various acoustic‐gravity wave modes due to volcano eruption were observed by accurate Beidou geostationary orbit (GEO) TEC, and the huge ionospheric hole was mainly caused by intense shock‐acoustic impulses. TEC rate of change index revealed globally propagating ionospheric disturbances at a prevailing Lamb‐wave mode of ∼315 m/s; the large‐scale EPBs could be seeded by acoustic‐gravity resonance and coupling to less‐damped Lamb waves, under a favorable condition of volcano‐induced enhancement of dusktime plasma upward E×B drift and postsunset rise of the equatorial ionospheric F‐layer. 

Abstract This study provides first storm time observations of the westward‐propagating medium‐scale traveling ionospheric disturbances (MSTIDs), particularly, associated with characteristic subauroral storm time features, storm‐enhanced density (SED), subauroral polarization stream (SAPS), and enhanced thermospheric westward winds over the continental US. In the four recent (2017–2019) geomagnetic storm cases examined in this study (i.e., 2018‐08‐25/26, 2017‐09‐07/08, 2017‐05‐27/28, and 2016‐02‐02/03 with minimum SYM‐H index −206, −146, −142, and −58 nT, respectively), MSTIDs were observed from dusk‐to‐midnight local times predominately during the intervals of interplanetary magnetic field (IMF) Bz stably southward. Multiple wavefronts of the TIDs were elongated NW‐SE, 2°–3° longitude apart, and southwestward propagated at a range of zonal phase speeds between 100 and 300 m/s. These TIDs initiated in the northeastern US and intensified or developed in the central US with either the coincident SED structure (especially the SED basis region) or concurrent small electron density patches adjacent to the SED. Observations also indicate coincident intense storm time electric fields associated with the magnetosphere–ionosphere–thermosphere coupling electrodynamics at subauroral latitudes (such as SAPS) as well as enhanced thermospheric westward winds. We speculate that these electric fields trigger plasma instability (with large growth rates) and MSTIDs. These electrified MSTIDs propagated westward along with the background westward ion flow which resulted from the disturbance westward wind dynamo and/or SAPS. 

Abstract We provide evidence that midlatitude postsunrise traveling ionospheric disturbances (TIDs) are comprised of electrified waves with an eastward propagation component. The post‐sunrise gravity wave (GW) wind‐induced dynamo action effectively generated periodic meridional polarization electric fields (PEFs), facilitating TID zonal propagation in a similar fashion as GW‐driven neutral perturbations. A combination of near‐simultaneous eastward and upward observations using the Millstone Hill incoherent scatter radar along with 2‐dimensional total electron content maps allowed resolution of TID vertical and horizontal propagation as well as zonal ion drifts(meridional PEFs). In multiple observations,oscillated in the early morning during periods when TIDs exhibited downward phase progression, 30–60 min period,140 m/s eastward speed, and 70 km vertical wavelength. Inside these TIDs, multiple flow vortexes occurred in a vertical‐zonal plane spanning the ionospheric topside and bottomside. Subsequently, PEFs weakened after a few hours as TID horizontal wavefronts rotated clockwise.